The event-plane method is obsolete 
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We argue that the traditional event-plane method, which is still widely used to analyze anisotropic 
flow in ultrarelativistic heavy-ion collisions, should be abandoned because flow fluctuations intro- 
duce an uncontrolled bias in the measurement. Instead, one should use an alternative, such as the 
scalar-product method or cumulant method, which always measures an unambiguous property of 
the underlying anisotropic flow and therefore eliminates this bias, and does so without any disad- 
vantages. It is known that this correction is important for precision comparisons of traditional v n 
measurements requiring better than a few percent accuracy. However, we show that it is absolutely 
essential for correlations between different harmonics, which can differ from the nominally-measured 
quantity by as much as 60% to a factor two, depending on the implementation, which is not standard 
between groups. We also describe how, using the corrected analysis method, the information from 
different subevents can be combined in order to optimize the precision of analyses. 



I. INTRODUCTION 

Anisotropic flow v n , in particular elliptic flow V2, is 
one of the most important observables of ultrarelativistic 
heavy-ion collisions p]-@. ^ is therefore important to 
devise the best possible methods to analyze it. A good 
measurement must be reproducible; in particular, it must 
be done in such a way that one can easily compare results 
from different experiments, using different detectors. For 
the sake of comparison with theory, it is also important to 
measure well-defined quantities, corresponding to generic 
properties of the system that are closely related to an 
interesting theoretical concept. 

In this paper, we argue that the event-plane 
method @, which is still used b y m ost of the major ex- 
perimental collaborations @, H, Tlfjl - fT^i does not satisfy 
any of the above requirements defining a good measure- 
ment. When the method was first developed, it was as- 
sumed that event-by-event fluctuations in v n were neg- 
ligible. However, it is now known that v n fluctuates 
significantly within a class of events fl3Hl5j . As a re- 
sult, an event-plane measurement of v n yields an am- 
biguous measure lying somewhere between the event- 
averaged mean value (v n ) and the root-mean-square 
value (v^) 1 / 2 [H, [13 ■ Where exactly depends on the 
"resolution", which strongly depends on the experimen- 
tal setup. This means that the same Pb-Pb collisions at 
the LHC yield a systematically different value of v n de- 
pending on whether they are analyzed by ALICE, CMS 
or ATLAS, and that comparisons to earlier measure- 
ments from lower energy collisions at RHIC are ambigu- 
ous. The difference is typically of a few percent for elliptic 
flow [HI , but is larger for higher harmonics like triangular 
flow i>3 [Hi, which is solely due to fluctuations. 

A slight variant of the event-plane method such as 
the scalar-product method introduced by the STAR col- 
laboratio n |18| or equivalently a two-particle cumulant 
method [13. |20|| make for a superior measurement be- 
cause they consistently yield the rms value of v n , while 
introducing no disadvantage compared to the traditional 
event-plane measurements. Supplemented by higher 



order cumulant measurements, these make traditional 
event-plane measurements redundant as well as ambigu- 
ous, as they contain no independent information about 
the underlying flow. 

Traditional event-plane analysis methods have also 
been used to measure higher-order, mixed-harmonic cor- 
relations |2lH23j |. These measurements also suffer from 
similar biases as traditional event-plane v n measure- 
ments, but in this case they are much larger 



24 , caus- 



ing significant problems even for qualitative comparisons. 
Fortunately, these problems can be repaired with a few 
simple changes. 

In Sec. UH we explain why the traditional event-plane 
method yields an ambiguous measurement. We show how 
the ambiguity can be simply removed. We also describe 
a new systematic way of improving the accuracy of the 
analysis by combining reference flows from "subevents" . 
In Sec. IIII1 the discussion is extended to the analysis of 
mixed correlations, where we demonstrate that the pro- 
posed corrected procedure is critical. 



II. MEASURING v n 

The picture underlying anisotropic flow measurements 
is that, if the system exhibits strong collective behavior 
(i.e., "flow"), particles in each event are emitted indepen- 
dently according to an underlying one-particle probabil- 
ity distribution [25| 
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where <fi is the azimuthal direction of an emitted particle, 
v n is the amplitude of anisotropic flow in the nth har- 
monic, and <£>„ the corresponding reference angle. Fur- 
ther assumptions underlying most analyses are that 
represents a global angle that depends little on pseudo- 
rapidity and transverse momentum, and that event-by- 
event fluctuations of v n also depend little on pseudora- 
pidity and transverse momentum. 
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Experiments typically detect anywhere from a few hun- 
dred particles to at most a few thousand in a given col- 
lision event, while the anisotropy coefficients v n are typ- 
ically a few percent or less. A precise reconstruction of 
the underlying probability distribution is therefore im- 
possible in a single event. Information about the under- 
lying probability distribution can only be extracted from 
(azimuthally symmetric) correlations between out goin g 
particles, averaged over a large ensemble of events |2q [. 

Historically, the most common way to access informa- 
tion about the coefficients v n has been so-called event- 
plane measurements. In these analyses, v n is extracted 
from correlations between particles of interest and a ref- 
erence detector which is typically separated by a gap in 
pseudorapidity [3, El]- In each event, one defines the flow 
vector Q n = (\Q n \ cos!//?*,,). |Q„j sin(??*„)) in harmonic 
n by @ 

\Q n \ cos(n\l> n ) = ^2 w j cos(n0j) 

3 

\Q n \ sin(n*„) = y^Wjsin(n0j) (2) 

3 

where the sum runs over particles seen in the reference 
detector, with azimuthal angles tf>j, and Wj is a weig htQ 
Alternatively, if the reference detector is a calorimeter, <f>j 
denotes the azimuthal angle of a segment of the calorime- 
ter, and Wj the energy deposited in the segment. 

The idea behind the event-plane method is that the 
direction of the flow vector \?„ provides an estimate of 
the corresponding angle <&„ in the underlying probabil- 
ity distribution [29( . Because a finite sample of particles 
is used, statistical fluctuations cause ^ n to differ from 
<!>„. The goal is to measure a quantity that is insensitive 
to these trivial statistical fluctuations and measures an 
unambiguous property of the underlying probability dis- 
tribution Eq. (JIJ. In the standard event-plane method, 
then, the final measurement is defined as 



v n {EP} 



(cosn(</>- 1> ra )) 
Res 



(3) 



where angular brackets denote an average over particles 
(with azimuthal angles <fi) and events. Res denotes a res- 
olution correction, which aims to correct for these trivial 
statistical fluctuations. To do this, one measures corre- 
lations between the reference detector in question and 
one or more additional reference detectors. There are 
essentially two cases, depending on the detector setup. 

f . If there are two identical reference detectors A and 
B symmetric around midrapidity and separated by 
a rapidity gap [|| (two-subevent method @, |29|). 
then 



Res ee v/(cosn(*„A - f„s)). 



(4) 



2. If there is no such symmetry, reconstructing the ref- 
erence flow requires information from a third [9| ref- 
erence detector C, where A, B and C are pairwise 
separated b y p seudorapidity gaps (three-subevent 
method 0, El GJ, HI) The resolution correction 
for v n measured with respect to the event plane 
from detector A is then 

Reg ^ j (cosn(V nA - ^ nB ))(cosn(y nA - V n ^)) 
CS ~ V (cosn(*„ B - *„c)) 

When the method was developed, it was assumed that 
dynamical fluctuations of the underlying probability dis- 
tribution are negligible. As such, the two underlying hy- 
potheses were: 

• The numerator of Eq. ((3]), (cosn(0 — $„)), factor- 
izes into (cosn((f>— $„)) (cos n($ n — ^ n ))- 

• The resolution correction Res is equal to 
(cosn($„ - *„)). 

In this case, it is clear that the resolution correction com- 
pletely cancels the statistical fluctuations coming from 
the use of a finite sample of particles in the reference 
detector, and v n {EP} measures directly (cos n(<f>— 3>™))- 

However, both hypotheses fail if v n in the underlying 
probability distribution fluctuates. The strength of the 
correlations between <j> and $„, and between and $„, 
both increase with v n , which causes factorization to break 
down. Similarly, the resolution correction is no longer a 
simple projection of the measured event plane ^> n onto 
the "true" event plane 

In the limit of perfect resolution (i.e., the number of 
particles used in the reference detector goes to infinity), 
^ n = $n a nd v n {EP} does indeed measure the event- 
averaged mean v n from Eq. ([TJ, 



v n {EP} -> (cosn(0- $„)) = (v n ). 



(6) 



In reality, the resolution is not perfect, and the result 
is usually closer to the low resolution limit [Hj]. In this 
limit, the numerator of Eq. is a two-particle correla- 
tion, which measures the average value of (v n v n ), where 
v n is the integrated v n in the reference detector. Sim- 
ilarly, the resolution correction is the square root of a 
two partic le correlation in the reference detector, that 
i s i V ■ The event-plane measurement thus yields a 
root-mean-square value, 



v n {EP} -+ Vra- 



(7) 



In general, the event-plane method yields a result 
which may lie anywhere between the two limits (|6]) and 



1 The simplest choices are Wj 
(Pt)/(Pt) forn = l E|. 



1 for n > 2, and Wj 



PT 



2 For simplicity we assume here that event-by-event flow fluctu- 
ations depend little on rapidity and transverse momentum, so 
that v n /v n is the same in every event. 
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FIG. 1. (Color online) Ratio of the rms v n , corresponding 
to the low-resolution limit of Eq. ((3} integrated over px, to 
the mean v n , corresponding to the high- resolution limit, as a 
function of the number of participant nucleons (binned in 5% 
centrality intervals), for a Pb-Pb collision at 2.76 TeV per nu- 
cleon pair. We model flow fluctuations by assuming v n oc e„ 
in each event, where e„ is the anisotropy of the initial distri- 
bution in the corresponding harmonic, defined as in Ref. [30| . 
The value of e n in each event is obtained from the Phobos 
Monte Carlo Glauber model [3l| . For 113, the ratio is close to 
~ 1.13, corresponding to Gaussian fluctuations [3^ ]. 



([7]). The exact measured quantity depends the detec- 
tor acceptance — even after the resolution correction is 
applied — so that the measurement is ambiguous. 

The difference between the two limits ranges from 6 
to 13%, as illustrated in Fig. [TJ As the field enters an 
era of precision phy sics, comparisons of different experi- 
mental analyses [17| and of theory to experiment [33| are 
significantly complicated by this dependence on analysis 
details. For example, the large detector acceptance and 
high multiplicities of recent measurements at the LHC 
allow for a larger resolution that was available at RHIC. 
The natural choice (and what has been done) is to take 
advantage of this extra resolution to minimize the statis- 
tical uncertainty of the measurement. The penalty, how- 
ever, is that comparisons to earlier measurements become 
ambiguous, and one has difficulty determining from these 
measurements how much the flow changes, e.g., with col- 
lision energy. Theoretical comparisons are hampered for 
the same reason. Theorists have access to more infor- 
mation than experiments, for example neutral particles 
in addition to charged. However, if they succumb to the 
temptation to make use of this extra information to im- 
prove statistics, comparisons to measurements become 
ambiguous and unreliable — even if they painstakingly 
follow the exact same procedures as the experiment [34j — 

m. 

Fortunately, this ambiguity can be removed by a sim- 
ple chan ge t o this procedure, known as the scalar-product 
method [l8j. Although it was originally introduced for 
unrelated reasons, it turns out to give a well-defined mea- 



surement in the presence of flow fluctuations [17| , despite 
differing only trivially from the above analysis. In this 
method, one simply replaces Eq. §5§ with 



v„{SP) = (IQ"l C0S "(^-^n)) 



Q, 



(8) 



where the resolution correction Q n is determined in the 
same way as the traditional event-plane method except 
without dividing by the magnitude of the measured flow 
vector \Q n \- 

Again there are essentially two cases. For identical and 
symmetric reference detectors we have 



QnA = V (QnA 1 QnB) 

while more generally one can use 



(9) 



r. _ (QnA ■ QnB) {QnA ' Qnc) , 1n , 

QnA=\ -=i =; ; ■ (10) 

V (QnB ■ Qnc) 

When flow fluctuations are negligible, the scalar prod- 
uct method thus defined is equivalent to the traditional 
event-plane method. When v n fluctuates, Q n A defined 
by either ([9]) or (fTUf scales like the rms value of v n , while 
the numerator of Eq. (jHJ) scales like (v n v n ). Therefore, 
the scalar-product method always yields the root-mean- 
square v n , regardless of the details of the analysis, and 
makes for a superior measurement. 

Equivalently, one can choose to measure a direct two- 
particle correlation [33] , which also always measures this 
same unambiguous root-mean-square value [1 3j j . It is im- 
portant to emphasize, however, that even when an event- 
plane- type analysis is preferred for practical reasons, only 
a trivial change to the standard analysis is required in or- 
der to obtain a well-defined observable. 

If the resolution is large enough, higher-order cumu- 
lants [l9|, H3] can be analyzed. These then yield an un- 
ambiguous measurement of higher-order, even moments 
of the distribution of v n , i.e. ((v n ) 2k ) [13j. Since a stan- 
dard event-plane method measures something different, 
one might naively think it useful to perform both analy- 
ses in order to obtain the maximum amount of indepen- 
dent information. However, it has been proven [ItJ that 
the event-plane method, in addition to being ambigu- 
ous, contains no independent information beyond that 
contained in the first two cumulants, v n {2} and i> n {4} 
(regardless of whether nonflow correlations are present 
in the system). So there is no reason to use a standard 
event-plane method in any future analysis. 

We end this section with a brief description of how 
the precision of the scalar-product method can be im- 
proved beyond previous implementations. Since at least 
two reference detectors A and B are needed in order to 
estimate the reference flow Qa, it is tempting to com- 
bine the information from both detectors into a single 
measurement with reduced statistical uncertainty. In the 
standard event-plane method, this can be done for identi- 
cal subevents at the expense of additional hypotheses and 
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algebraic complications [Qj, i.3f§| , but there is no known way 
of recombining non-identical subevents, as in the three- 
subevent method. 

With the scalar-product method, combining reference 
flows from two reference detectors A and B is straight- 
forward, even if they are not identical. One simply mea- 
sures v n independently with respect to A and to B. By 
linearity of the scalar product, combining both measure- 
ments (denoted by v„,a and i> n ,s) amounts to taking a 
weighted average. We show in Appendix [A] that the op- 
timal weighting is 



XA v n,A + XB V n,B 

Xa + X 2 b 



(11) 



where \A is the dimensionless resolution parameter [9], 



XA 



InA 



((Qn,A) 2 ) ~ QIa 



(12) 



When detectors are combined, their resolution parame- 
ters x a dd in quadrature. Generalization to three refer- 
ence detectors A, B, C is straightforward. 



III. HIGHER-ORDER CORRELATIONS 

Fourier harmonics of the azimuthal distribution are 
now measured up to i>4 [l(| at RHIC and up to v$ @, Ho| 
and vq 7} at LHC. Much additional information is con- 
tained in mixed correlations between different harmon- 
ics [4H - I44I ] , which give information about event plane an- 
gles and will likely play a large part in upcoming 
measurements. Before rushing into analysis, it is impor- 
tant to define observables for good measurements, which 
will allow unambiguous comparison between different ex- 
periments, and with theory. 

An early measurement of this type, performed before 
the importance of flow fluctuations was noticed, is V4 with 
respect to the event-plane of elliptic flow *§> 2 [SI, [22j]. 
The high resolution limit of these measurements is the 
nominal quantity 



v 4 {^ 2 } (V4 cos 4($ 4 - $ 2 ))- 



(13) 



The actual measured quantity again depends on details 
such as the event-plane resolution. In addition, there is 
an ambiguity coming from differences in the implementa- 
tion of the "resolution correction" . The seminal paper Q 
proposes two possible implementations, which are equiv- 
alent when there are no flow fluctuations, but turn out 
to be different in the presence of fluctuations. Thus one 
can obtain either of two expressions for the low- resolution 
limit of the measurement, depending on the actual im- 
plementation. STAR HH has used Eq. (11) of Ref. d, 
which reduces to 



v 4 {^ 2 } 



(vivl cos4($ 4 - $2)) 

W) 



(14) 



in the low-resolution limit, where again v 2 is integrated 
over the reference detector. PHENIX instead has used 
Eq. (14) or Eq. (16) of 9], which yields, in the low- 
resolution limit, 



V 4 {^ 2 } 



(u 4 uf cos4($ 4 - $2)) 

7W) ' 



(15) 



The actual measurement should lie between the high- and 
low-resolution limits corresponding to the respective im- 
plementation. As recently demonstrated in Ref. [33[ in 
full event- by- event hydrodynamic calculations, the de- 
pendence on resolution is larger than for the v n {EP} 
measurements described above. In fact, this should al- 
ready have been apparent. In addition to v 4 {ty 2 }, the 
STAR collaboration also measured a 3-particle correla- 
tion t>4{3} [2lj . which corresponds unambiguously to the 
low-resolution limit of their event-plane measurement. 
They found that U4{3} is larger than ^{vf^} by 20% 
for mid central collisions. This difference is much larger 
than the difference between, e.g., v 2 {EP} and the low 
resolution limit v 2 {SP} in the same experiment (which 
is roughly 5%). 

Further, the difference in implementation, coupled 
with flow fluctuations, is actually the most likely ex- 
planation for the puzzling fact that the STAR measure- 
ments are systematically larger than the PHENIX mea- 
surements of what was claimed to be the same quan- 
tity [22|, [24]]. The fact that this may not have been 
fully understood until now underlines the problems aris- 
ing from ambiguities in experimental analyses combined 
with the presence of flow fluctuations. 

Recently, the ATLAS collaboration has released pre- 
liminary measurements of an extensive set of correlations 
between event planes [23[ in Pb-Pb collisions at the LHC, 
for which they decided to use traditional event-plane 
methods. To illustrate the large ambiguities present in 
these measurements, we consider the first of these corre- 
lations, which is similar to the earlier measurements of 
V4{^> 2 }. It is a correlation between the event planes ^4 
and ^2, namely, (cos4(\E r 4 — ^ 2 )), which is then divided 
by a resolution factor that is supposed to to unravel the 
correlation between the underlying reference directions 
$4 and $2- They thus write: 



(cos4($ 4 - $2)) = 



(cos4(* 4 - vJ/ 2 )) 



(cos4(* 4 - $4)) (cos 4(^2 - $2)) ' 

(16) 

However, this equation is based on the same hypothe- 
ses as the event-plane method, which fail if v 2 and V4 
fluctuate. Specifically, the numerator of the right-hand 
side does not factorize into (cos 4(^4 — $4)) (cos 4(\E , 2 — 
<I>2))(cos4($4 — $2)), and the resolution corrections esti- 
mated using standard procedures do not actually corre- 
spond to (cos4(* 4 - $ 4 )) and (cos4(* 2 - $a))< 

Specifically, the resolution correction (cos 4(^4 — $4)) 
is defined by Eq. (j4|) or ([5]). The resolution correc- 
tion (cos4('I'2 — $2)) corresponds to the convention that 
PHENX used for 1*4 {"J^}, which is defined by a simple 
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generalization of Eq. (HJ (or ([5])), namely, 

(C0S4(# 2 - $2)) = yJ{cOS^2A ~ *2S))- (17) 

In the limit where the resolutions on both harmonics 
are perfect, the quantity measured by ATLAS is indeed 
the nominal value, corresponding to the left-hand side of 
Eq. d]): 

cos4(* 4 - ^2){EP} ->■ (cos4($ 4 - $ 2 ))- (18) 

In the limit of low resolution, however, repeating the 
same reasoning as in the previous section, what is ac- 
tually measured is the quantity 

00B4(* 4 -* a){ EP}-,M^^M. ( 19 ) 

V( V i)(«2> 

As before, the measured value is generally between the 
two limits0 Similar expressions can be written for the 
other correlations in the large set measured by ATLAS. 
Unlike for measurements of v n , however, there is no 
known way to even estimate where between these lim- 
its the result lies based only on the reported event-plane 
resolution. 

In order to illustrate the dependence of these measure- 
ments on resolution, we use here a simplified calculation, 
though the results are consistent with more sophisticated 
calculations (33j . It is known from event-by-event hydro- 
dynamic calculations that pr-integrated t>2, v 4 , and the 
corresponding reference directions $2 and $4 can be ac- 
curately predicted in each event from the initial eccen- 
tricities according to [3(| 

v 4 e il ^ = 6e 4 e 4i *'"' 4 + c(e 2 e 21 *"- 2 ) 2 , (20) 

where a, b, c are real coefficients, and e„ and Qi n ,n are 
the initial anisotropy in harmonic n and its reference di- 
rection, as defined in Ref. (30| . 

We thus compute e n and §i n ,n using the PHOBOS 
Monte Carlo Glauber [3l| and take the values of b and 
c from the viscous hydrodynamic calculation by Teaney 
and Yan [3| , for simplicity taking centrality-independent 
values of b = 0.018 and c = 0.070. 

The right-hand side of Eqs. (fT8|) and (fT9|) are displayed 
in Fig. [5] (top, left), as well with their ratio (bottom, 
left), and similarly for Eqs. ((131), (JH) and (TT5) on the 
right side of the figure. The ratio varies from 1.2 to 2 (in 



3 Note that the low-resolution limit 1 1191) lies between —1 and +1 by 
construction, even though it does not reduce to a simple angular 
correlation. 

4 Teaney and Yan use cumulants of the initial distribution, which 
are linear combinations of our moments: with their notation, 
b = w 4 /Ci and c = 104(22) /^l + 3{w 4 /C 4 )(r 2 ) 2 / (r 4 ) . We take 
U14/C4 ~ 0.018 and w 4 ( 2 2)/ £ 2 — 0-04 from their Fig. 2, right, 
and we estimate (r 2 ) 2 / (r 4 ) ~ 0.6 from a Glauber calculation. 
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FIG. 2. (Color online) Top, left: low and high resolution lim- 
its of the correlation between $2 and $4, defined by Eqs. fj 19p 
and QS) , versus number of participants. Top, right: low 
and high resolution limits for integrated ^{^j, defined by 
Eqs. (|13[) . (|14[) and (|15[) . Bottom panels: ratios of the low- 
resolution to high-resolution results. The nominally-measured 
values correspond to the high-resolution limit, which can dif- 
fer from what is actually measured (usually closer to the low- 
resolution limit) by up to a factor 2. All results are for Au+Au 
collisions at top RHIC energy. 

agreement with full event-by-event hydrodynamic calcu- 
lations of U4{^2} [iH ). Comparing with Fig.[TJ one sees 
that the ambiguity induced by the event-plane method 
is much worse for mixed correlations that for individual 
measurements of v n . 

As a recent illustration of the problems that can arise 
from this, we point to the correlations calculated in 
Ref. [IK, corresponding to the high resolution limit of 
the recent preliminary ATLAS measurements — i.e., 
the quantities that were nominally claimed as measured. 
They are all clearly smaller in magnitude than the mea- 
surements, despite each observable having the correct 
sign and centrality dependence. Naively, this seems to 
indicate that the theory needs to be changed somehow 
to accommodate this difference, but in fact the finite res- 
olution of the experiment could in principle explain the 
entire discrepancy. Further, while it is straightforward to 
calculate either the low or the high resolution limit, as 
done in that work, a calculation that reliably corresponds 
exactly to the measured quantity adds considerable dif- 
ficulty. In fact, it is not possible without access to more 
information about the experimental details than is avail- 
able. 

Fortunately, the ambiguity of the event-plane method 
can easily be removed by using a straightforward gen- 
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eralization of the scalar-product method. We write the 
flow vector @ in complex form, which allows to write 
compact expressions for mixed correlations: 



(21) 



With this notation, one can unambiguously obtain the 
right-hand side of Eq. (fl"9]l by defining 



(«4«2 cos4($ 4 - $ 2 )) 



Re((Q 2 A) 2 Ql B ) 



^fWM ' ^(q±aQi b )((Q2aQ* 2B ) 2 Y 

(22) 

where Re stands for the real part, subscripts A and B 
refer to two symmetric reference detector separated by a 
rapidity gap, and Q* nB denotes the complex conjugate of 
Q n B- In the same way as for individual measurements 
of v n , this method yields a result which is unambiguous 
and independent of the number of particles selected for 
analysis. 

Generalization of the above discussion to other mixed 
correlations [H, EJ is straightforward. Specifically, the 
ATLAS collaboration measures correlations of the type 
(cos(ci5'i + • • • + Ici^i)), which are divided by the prod- 
uct of resolution factors in each harmonic (Eq. (1) of 
Ref. [23j]). Ambiguity is removed by inserting appropri- 
ate factors of \Q n \. That is, the numerator should be 
replaced as 



(cos(ci^i 

-> (\Qi 



+ lci^i)) 

|Q i |l C! lcos(c 1 * 1 + --- + Z Q ^)), (23) 



while the resolution factor in harmonic j, denoted by 
Res{cjj^J j}, should be replaced with a straightforward 
generalization of Eq. (ITU|) : 



^){{QiAQ* jC r A ) 



((QjbQ* c )^) 



(24) 



In the case of two identical and symmetric reference de- 
tectors, this can again be simplified to use 



Res{ Cj j$ ,} 



((Q. 



(25) 



With this convention for the resolution factors, the 
Cauchy-Schwarz theorem ensures that all the two-plane 
correlations lie between —1 and +1, though there is no 
such restriction for three-plane correlations (such as the 
correlation between $2, $3 and $5) which could in prin- 
ciple exceed unity. Another option is to scale the numer- 
ator as proposed in Ref. 41], which amounts to scaling 
Eq. ([23]) by Qj 01 ' • • • Q^ J ' , where the resolution factors Qj 
are given by Eq. (j9|) or (flO|) . All possible choices are re- 
lated to each other by information contained in standard 
cumulant measurements of v n , so it is only necessary to 
choose one convention. The only important considera- 
tion is to choose something that does not contain the 
ambiguities of a traditional event-plane analysis. 



Note that the particular correlation between $2 and 
$4 studied above can be contaminated by nonfiow corre- 
lations; Therefore it is important that the flow vectors Q2 
and Q4 are determined in regions separated by a pseu- 
dorapidity gap 23]. Other mixed correlations, such as 
the correlation between $2 and $3, are not sensitive to 
nonfiow correlations (though the resolution factor in the 
denominator still requires subevents that are separated 
in pseudorapidity). The need for gaps in pseudorapidity 
must be considered on a case-by-case basis [46[ . 



IV. CONCLUSION 

We have explained that measurements of anisotropic 
flow using the event-plane method are ambiguous be- 
cause the result depends systematically on analysis de- 
tails. This dependence, which is due to event- by- event 
flow fluctuations, is typically up to a few percent for v-i 
and 10% for U3 and higher harmonics. For mixed corre- 
lations involving event planes from different harmonics, 
the dependence is much larger — we have demonstrated 
that the result may vary by as much as a factor 2 de- 
pending on analysis details. We have shown that these 
ambiguities of the event-plane method can be easily re- 
paired with only minor modifications and without any 
additional complication or drawbacks. We therefore con- 
clude that the traditional event-plane method should be 
abandoned in all future analyses. 
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Appendix A: Combining reference flows 

In this Appendix, we show that Eq. (ITT|) defines the 
optimal combination between anisotropic flows measured 
in detectors A and B using the scalar-product method, 
in the sense that it minimizes the statistical error. 

Eq. ([5]) can be rewritten as 



V n {SP} 



(Qn ■ U) 
Qn 



(Al) 



where u = (cos ruj>, sin ncj)) is the unit vector of the par- 
ticle in harmonic n. Assuming that v n -C 1, so that u is 
essentially a random unit vector, and assuming that one 
measures v n in a narrow bin which contains at most one 
particle per event, a straightforward calculation shows 
that the statistical error on v n {SP} is (for simplicity, we 



7 



drop the index n from now on) 



Sv 



1 



(Q 2 ) 



/2N Q 



(A2) 



where N is the number of particles in the bin. 

Consider two subevents A and B, separated by a gap 
in pseudorapidity. We denote by Qa and Qb the corre- 
sponding flow vectors in harmonic n, and by Qa and Qb 
the corresponding resolution corrections. If one measures 
v n in a detector that is separated in pseudorapidity from 
both A and B, the measurement can be done using either 
A or B as a reference: 



va.b 



(Q 



A,B ■ U) 



Q 



A.B 



(A3) 



One can combine linearly the flow vectors and the reso- 
lution corrections 



Q = Q A + XQb 
Q = Qa + XQb, 



(A4) 



where A is arbitrary. The resulting combined measure- 
ment is 



^combined 



_ (Q ■ u) Qava + XQbvb 



Q 



Qa + XQb 



(A5) 



Using Eq. (IA2[) . straightforward algebra shows that the 
value of A which minimizes the statistical error on 

^combined IS 



A = 



Qb(Qa) " Qa(Qa ■ Qb) 
Qa(Q 2 b)~Qb(Qa-Qb) 



(A6) 



Since (Qa ■ Qb) = QaQb by definition, this can be 
rewritten as 



Qb((Q 2 a)-Q 2 a ) _ QaxI 



(A7) 



where the resolution parameter xa has been defined in 
Eq. (fT"2|) . Inserting this value of A into Eq. (IA5|) . one 
recovers Eq. (TTT1) . 
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